Air-clad fibers: pump absorption 
assisted by chaotic wave dynamics? 

Niels Asger Mortensen 

MIC — Department of Micro and Nanotechnology, Nano»DTU, Technical University of 
Denmark, 0rsteds Plads, DTU-building 345 east, DK-2800 Kongens Lyngby, Denmark. 

nam@mic. dtu. dk 
http://www.mic.dtu.dk/nam 



Abstract: Wave chaos is a concept which has already proved its 
practical usefulness in design of double-clad fibers for cladding-pumped 
fiber lasers and fiber amplifiers. In general, classically chaotic geometries 
will favor strong pump absorption and we address the extent of chaotic 
wave dynamics in typical air-clad geometries. While air-clad structures 
supporting sup-wavelength convex air-glass interfaces (viewed from the 
high-index side) will promote chaotic dynamics we find guidance of regular 
whispering-gallery modes in air-clad structures resembling an overall 
cylindrical symmetry. Highly symmetric air-clad structures may thus 
suppress the pump-absorption efficiency 77 below the ergodic scaling law 
77 <*: A c /A c \, where A c and A c \ are the areas of the rare-earth doped core and 
the cladding, respectively. 
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1. Introduction 

Cladding-pumped fibers are in general considered interesting candidates for use as high-power 
fiber lasers and amplifiers JTJ [2j [3] [4J . Fiber lasers provide the advantages of high power levels 
combined with very high beam quality. The successful scaling of power levels for fiber lasers 
is powered by an ongoing technological progress along several lines including the development 
of diode lasers, pump coupling schemes, and improved double-clad fibers. 

The inset in Fig. [T] shows a schematic presentation of a cladding-pumped fiber laser. The 
double-clad approach was pioneered by Snitzer and co-workers in the late 1980's (see Ref. J5) 
and references therein) and the basic concept of the double-clad fiber laser is to facilitate con- 
version of multi-mode light into a single-mode beam of high quality. The pump light is corralled 
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Fig. 1. Plot of the pump-absorption efficiency T](L). The solid line shows the ergodic limit, 
Eq. (0, while the dashed line indicates the initial dynamics and the critical length L*, 
Eq. i ll It . The inset illustrates a setup employing a double-clad fiber for conversion of multi- 
mode light with a low beam quality into single-mode light with a high beam quality. 



Fig. 2. Panel (a) shows a micrograph of an air-clad photonic crystal fiber made in silica 
(courtesy Crystal Fibre A/S). The pump light is corralled by an air clad consisting of 42 
air holes (dark regions) separated by sub-micron silica bridges. The diameter of the inner 
cladding is of the order 125 flm. Panels (b) to (h) show time-dependent finite-element sim- 
ulations of a scalar wave equation at different times (click panel to view animation, 3.3 
Mbyte) for the initial Gaussian excitation shown in panel (b). 

by the double-clad and while light propagates down the fiber, photons are absorbed and con- 
verted to stimulated emission by a rare-earth doped fiber core at the center of the inner cladding. 
Typically, the laser cavity is formed by inscribing Bragg reflectors in the fiber core. 

Here, we focus on air-clad technology H which has greatly advanced the potential of double- 
clad fibers Q |U |9] [TO] QT| [T2] . The large index contrast between air and silica makes the 
cladding support pump light of very high numerical aperture 1131 . thus facilitating pumping 
by relatively cheap high-NA pump sources. For the corralled pump light, skew rays lfT4l . or 
whispering gallery modes (WGMs), are well-known obstacles to efficient pump absorption 
since they do not have a significant overlap with the core. Such modes are particular inherent 
to cylinder symmetric fiber geometries, such as standard MCVD-fabricated fibers. The prob- 
lem may be circumvented by breaking the high cylindrical symmetry and D-shaped cladding 
geometries (T51 have become a standard approach to suppression of skew rays. 

Pump absorption in air-clad fibers is typically very high and it has been speculated [ 16 1 that 
this is coursed by chaotic ray dynamics. Non-integrable air-clad geometries seem to be inherent 
to the fabrication method where the shapes of the air holes arise in an interplay and competition 
of e.g. the glass viscosity and glass-air surface tension. Fig. [2] illustrates one example of an 
air-clad fiber iflOl which has a classically chaotic geometry with a convex boundary (from the 
point of view of the pump light) which will scatter incident waves in a highly non-deterministic 
way. Time-dependent finite-element simulations of a scalar wave equation in a corresponding 
confined geometry, panels (b) to (h), do indeed support this picture of chaotic wave dynamics. 
Indeed, the pump absorption observed in Ref. IflOl is much stronger than for an equivalent 
circular geometry. 

In this work we use finite-element simulations to show how the wave dynamics has finger- 
prints inherited from the classically chaotic geometry. In particular, we consider the statistics 
of the mode-overlap with the core and argue that air-clad fibers tend to be chaotic by nature in 
the case of sup-wavelength convex air-clad features. 



2. Formalism 



For simplicity we consider the simplified problem with a spatially homogeneous inner cladding 
region £l c \ with dielectric function e c i = and a rare-earth doped fiber core region i2 c with 
dielectric function £ c , see panel (a) in Pig. EJ The inner cladding may be defined by several air 
regions il a i, -. The complex- valued dielectric function is then of the form 



e(r) = e'(r) + ie"(r) = < 



, re£l cU (1) 
1 , r e Oajr, 



where for the rare-earth doped core e" <C e' c so that absorption is a weak perturbation. In the 
following we consider temporal harmonic modes which are governed by the wave equation 

VxVx|£ m )= e (^) 2 |£ m ), m = 1,2,3,... (2) 

where co„, is the angular frequency and c is the speed of light. For a fiber geometry, translational 
invariance along the fiber axis results in modes of the plane-wave form e ! (^~ ffl ') with the disper- 
sion relation co(/3 ) found by solving Eq. (|2j. Typically, the dispersion properties are obtained 
through numerical solutions of Eq. ©, but perturbation theory is another highly valuable tool 
in analyzing the consequences of minor perturbations such as variations Ae in the dielectric 
function. For a fixed frequency, we apply standard first-order perturbation theory to predict the 
corresponding shift Aj3 in the propagation constant j3, 

. _ / dtO m V 1 co m (E m \Ae\E m ) 



dp J 2 (E m \e\E m ) 

where v g = dco/df} is the group velocity which is mathematically introduced through the use 
of the chain rule ifTTl . 

3. Pump absorption 

Obviously, a complex dielectric function as in Eq. (Q~|) will for a fixed and real-valued frequency 
result in a complex valued propagation constant j5 = /3' + /J3" . The interpretation of the imag- 
inary part as a damping becomes obvious when expecting the intensity I(z,t) « \e'(P z ~ m ^ | 2 = 
e~ 2 ^ z = e~ az where we have introduced the damping parameter a = 2)3". In the following we 
will consider the case where the imaginary part in Eq. ([!} is a small perturbation to the real 
part, i.e. e" <C e' so that j3" <C j3'. From Eq. (O we immediately get 

co e" 

a m = 2ft,', =f m x j ~ f m x a c (4) 

Vg,m £ c 

where we have introduced the relative optical overlap with the core 

_ (E m \E\E m ^ ^ (E m \D m } c ^ 
(E„,\e\E m ) c+d (E m \D m ) c+cl 

with |Z)) = £\E} being the displacement field. In the last equality of Eq. we have approxi- 
mated the group velocity by that in the homogenous core material, thus allowing for the intro- 
duction of the core material absorption parameter a c . The interpretation of Eq. is straight 



forward; the absorption of a given mode is very intuitively given by a c of the core material 
weighted by the relative optical overlap / of the mode with the core. Now, since the cladding 
diameter by far exceeds the wavelength of the light we will have a continuum of cladding 
states, even for typical low-NA air-clad fibers ifTSl . so studying the characteristics of the indi- 
vidual modes is cumbersome. Instead it is common to quantify pump absorption by a single 
parameter, i.e. the pump absorption efficiency as indicated in Fig.Q] 

n{L)=f daP(a)[l-exp(-aL)}~ f dfP{f) [1 -exp(-/ x cfcL)] . (6) 
Jo Jo 

Here, P(a) is the distribution of attenuation coefficients while [1 — exp (— uL)] is the absorption 
efficiency for the individual mode for a length L of propagation. In the second equality we have 
used the approximation in Eq. <j4j and introduced the distribution P(f) of the overlap integrals 
defined by Eq. (|5). With the form of Eq. (O we have implicitly assumed a uniform excitation at 
the input facet of the fiber so that each mode has the same initial amplitude and initially carries 
the same power. Furthermore, mode mixing due to longitudinal non-uniformities, including 
bending, is obviously not taken into account. In that sense Eq. © gives a lower bound for the 
pump absorption efficiency. 

A general evaluation of the integral of course requires that we know the distribution func- 
tion in detail. However, the initial decay is easily found by expanding the square-bracket 
part of the integrand. This conveniently leaves us with a result in terms of the moments 
(/") = Jo ' d fP{f)f n of the distribution function, 

tj(L) - (/) a c L- 1 - (f) (a c L) 2 + ^([a c L] 3 ). (7) 

Obviously, the higher a mean value (/} and the smaller a spread (/ 2 ) the more efficient a pump 
absorption for a given length L. From Eq. (O we may now easily estimate the critical length of 
fiber L* required to achieve an efficient conversion of the pump. From rj (L* ) ~ 100% we get 

L*~(f)- l xa- 1 (8) 

demonstrating how fiber geometries with a large mean value (/) will facilitate fiber lasers with 
a shorter laser cavity. In the following we will study chaotic means for enhancing the mean 
value and thus decreasing L*. 

4. Chaos and the ergodic hypothesis 

The term 'chaos' is used in many contexts of mathematics and physics, but here it is used with 
reference to the wave dynamics in a confined geometry. Since the double clad structures that 
we are considering have dimensions much larger than the wavelength it is quite natural to illus- 
trate the signatures of chaos by means of the classical ray dynamics that occurs if we consider 
geometrical optics. In Fig.[3]we illustrate the important difference between regular integrable 
geometries and their non-integrable chaotic counterparts. The circular geometry in panel (a) is 
an excellent example where (assuming total internal reflection) we may easily advance the rays 
in time, keeping in mind that the angle of reflection equals the angle of incidence. Panel (b) 
shows the corresponding whispering-gallery like ray pattern that results at a much later time. 
From a pump-absorption point of view the problem is of course that the rays do not visit the 
center of the geometry (where the absorbing core is typically located) and changing the initial 
conditions (by e.g. changing the angle of incidence) does not necessarily change this picture 
significantly. Panel (c) shows a non-integrable chaotic geometry where the inclusion of small 
circular scatterer completely changes the ray dynamics compared to panel (b) even though the 



Fig. 3. Ray-tracing dynamics illustrating the different response of integrable geometries, 
panels (a,b), and non-integrable chaotic geometries, panels (c,d). Panels (a) shows initial 
ray trajectories in the transverse section of an integrable geometry corresponding to an 
idealized multi-mode fibers having perfect circular symmetry. Panel (b) shows the regular 
whispering-gallery like pattern that results after propagating the initial rays in panel (a) for 
a long time (> 10 3 round trips). Panel (c) shows corresponding results for a non-integrable 
chaotic geometry including a circular scatterer. Note that the ray pattern results from the 
same initial conditions as for panel (b). Panel (d) shows a corresponding 'scarred' wave 
function obtained from a finite-element solution of a scalar wave equation. 



two trajectory patterns originate from the same initial conditions. If we focus on the wave na- 
ture of light in classically chaotic geometries then chaos reveals itself in several ways including 
the 'scarred' nature of the wave functions |fl9l as illustrated by the example in panel (d). As 
for the ray dynamics in panel (c) the wave functions also tend to explore the entire phase space 
rather than being of the whispering-gallery like nature seen in panel (b). 

From a pump absorption point of view, chaos completely changes the distribution P(f) 
of overlaps. Assuming that the eigenfields are, on average, spatially fully homogenized, then 
Eq. (0 reduces to 

where A c = dr and A c \ = dr are the areas of the core and cladding regions, respectively. 
In the following we use the definition fo = A c /A c \. The above assumption corresponds to the 
ergodic limit where all modes visit all parts of the phase space equally. We then have P(f) ~ 
8(f — fo) so that (f \ ~ fo and due to the simple properties of the Dirac delta function Eq. (O 
is now easily evaluated with the result 

r](L) ~ 1 - exp(-/oa c X) (ergodic limit) (10) 

showing a single-exponential behavior as illustrated in Fig.Q] Finally, Eq. ([H} reduces to 

L -/ l xa c 1 --xa c 1 (11) 

Ac 

which is the ultimate scaling which can be achieved for the absorption of the pump. 



5. Overlap statistics 

In the following we illustrate the difference between an integrable geometry and a chaotic 
geometry by comparing rods of glass in air with a circular and D-shaped cross section, respec- 
tively. We use a finite-element method (COMSOL Multiphysics) to calculate fully vectorial 
eigenmodes of Eq. (|2} from which we evaluate Eq. (0 numerically. Grouping the data into bins 



3 

2.5 

2 

S 1-5 

a, 

l 

0.5 



12 3 4 

f/fo 

Fig. 4. Distribution functions for a regular integrable and non-integrable chaotic geometry. 
The dashed lines shows numerical fits to Eq. H2i with (/) //o ~ 9% and (/)//o ~ 100% 
for the circular and D-shape geometries, respectively. 

of width /o/4we arrive at the histogram shown in Fig. |4] For simplicity we have considered 
£ c i = £ c = (1.45) 2 and for the wavelength we have used X/R =0.1 with R being the radius 
of the cladding. Furthermore, the core is placed in the center of the fiber and A c /A c \ = 0.01. 
Results for the circular case are shown in blue while the red results are for a corresponding D- 
shaped geometry. The dashed lines are numerical fits to the following normalized distribution 
functions 

f </> _1 «p (-//</» , </>//o~9% , (circular), 
Puif) = < (12) 
[ 4/(/)- 2 exp(-2//(/)) , </)//„ -100% , (D-shape), 

with the mean value (/) used as fitting parameter. Fig.|4]is illustrating the close to exponential 
distribution in the case of the integrable geometry while for the chaotic geometry the weight is 
shifted to much larger values with a suppression of modes with a poor overlap. The insets show 
typical modes for the two geometries. For the circular geometry the peak in the distribution 
originates from WPGs which are suppressed in the D-shaped geometry where modes are of the 
'scarred' nature with an average overlap (/) — /q in agreement with the ergodic hypothesis. 
Obviously, in the context of 77 the chaotic structure is superior to the integrable one as also 
emphasized previously in the literature iflBl l20l [2T| l22l . 

A quite natural question now arises; are air clad structures truly chaotic? As we shall see the 
answer to this inquiry depends a lot on the geometry into which the air holes are arranged. For 
the simulations in Fig. [2] we neglected the sub-micron silica bridges separating the sup-micron 
convex air-clad features. At first this appears reasonable if the bridges are sub-wavelength. 
However, one should not compare to the free-space wavelength A (or X /n c \ for that sake), 
but one should rather consider the effective transverse wavelength X± which may exceed A 
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Fig. 5. Qualitative illustration of the distribution function for an air-clad structure display- 
ing partly integrable (A^ 3> 8) and partly chaotic dynamics (A^ <C 5). The insets illustrate 
the two classes of modes representing a vanishing overlap / ~ and a finite overlap / ~ /o, 
respectively. 

significantly. Projecting the problem onto the transverse cross section of the fiber we get 



where n e ff = cj3 /©,„ is the effective index of mode m. Now, let 8 be a characteristic length scale 
for the air-clad features, such as the bridge width or the air-hole diameter. The modes may now 
be classified by comparing A^ with 8; 

• For X±^> 8 the electromagnetic field effectively experiences a smooth boundary defined 
by the air clad and the wave dynamics may be integrable or chaotic depending on the 
effective shape of the corral. 

• For Ax <C 8 the field explores all the detailed features and in particular the convex air-clad 
boundary so that the wave dynamics turn chaotic. 

While the geometry in Fig. [2] is chaotic in both limits one can easily imagine more regular 
air-clad structures which will only be 'partly' chaotic, i.e. low-NA modes will experience an 
integrable confinement geometry while high-NA modes will experience a chaotic geometry. 
Indeed, the insets in Fig. illustrate this perfectly for a geometry very much resembling the 
over-all cylindrical symmetry of the air-clad structure employed in more recent works such 
as Ref. |23|. The air-clad structure consists of 27 air-holes distributed evenly in a cylindrical 
symmetry with radius R. The modes are obtained by fully-vectorial finite-element solutions of 
Eq. (|2]i for a wavelength of X/R = 0.01. While a quantitative numerical simulation of the full 
statistics and the distribution function is cumbersome we imagine that the two classes of modes 
will give rise to a distribution function qualitatively interpolating the two limits in Eq. (fT2l 
as indicated schematically in Figj5] The existence of skew rays and WPGs, represented by 
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Fig. 6. Panel (a) shows a whispering-gallery mode in a ring-shaped air-clad geometry while 
panel (b) shows a corresponding mode for a D-shaped air-clad geometry. 

the exponential peak near / ~ 0, will eventually suppress the average value (/) below the 
ergodic limit given by /o thus compromising the pump-absorption efficiency 77 and the critical 
absorption length L* . Deforming the ring-shaped air-clad slightly, or introducing an additional 
air hole as in panel (c) or (d) of Fig. [3] would suppress the exponential peak at / ~ and 
restore the distribution function with (f \ ~ f characteristic of chaotic wave dynamics. Fig. [6] 
illustrates how WPGs may be suppressed by constructing the air-clad in a D-shaped fashion. 

6. Conclusion 

We have discussed wave chaos in the context of the pump-absorption efficiency 77 in air-clad 
fibers which are now being widely used for high-power cladding-pumped fiber lasers. Starting 
from the wave equation we rigorously derive a result for 77 which in the ergodic limit scales 
as 77 oc A c /A c \, where A c and A c \ are the areas of the rare-earth doped core and the cladding, 
respectively. While air-clad structures supporting sup-wavelength convex air-clad features will 
promote chaotic dynamics we also show how this may be jeopardized by whispering-gallery 
modes in air-clad structures resembling an overall cylindrical symmetry. In the absence of 
mode-mixing mechanisms, highly symmetric air-clad structures may thus suppress the pump- 
absorption efficiency 77 below the ergodic scaling law. 
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